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ABSTRACT 

The stationary hydrodynamic equations for transonic viscous accretion discs in Kerr 
geometry are derived. The consistent formulation is given for the viscous angular mo- 
mentum transport and the boundary conditions on the horizon of a central black 
hole. An expression for the thickness of the disc is obtained from the vertical Euler 
equation for general accretion flows with vanishing vertical velocity. Different solu- 
tion topologies are identified, characterized by a sonic transition close to or far from 
the marginally stable orbit. A numerical method is presented that allows to integrate 
the structure equations of transonic accretion flows. Global polytropic solutions for 
the disc structure are calculated, covering each topology and a wide range of physi- 
cal conditions. These solutions generally possess a sub-Keplerian angular momentum 
distribution and have maximum temperatures in the range fO 11 — fO 12 K. Accretion 
discs around rotating black holes are hotter and deposit less angular momentum on 
the central object than accretion discs around Schwarzschild black holes. 
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1 INTRODUCTION 

The presence of hot optically thin gas in accreting black hole systems such as X-ray binaries or active galactic nuclei has been 
established by X-ray and 7-ray observations. The spectra are very hard, with power-law photon indices on ~ 1.5 — 2.5 in the 
energy range E ~ keV to 100 keV (Gilfanov et al. 1995; Maisack et al. 1993; Johnson et al. 1994). They are usually fitted 
by a Sunyaev-Titarchuk (198C) thermal comptonization model with optical depth r cs <; 1 and electron temperature T c ~ 10 9 
K. These hard spectral components exclude as an interpretation all optically thick accretion models, like the standard thin 
accret ion disc model of Novikov & Thorne (1973; hereafter NT73) or the 'slim accretion disc' model of Abramowicz et al. 
( 1988 ), both of which cannot attain temperatures higher than a few ~ 10 7 K. The first hot optically thin disc model in 
agreem ent with hard spectral components was the two-temperature accretion disc model by Shapiro, Lightman & Eardley 
( 1976 ), with electron temperature T c ~ 10 8 — 10 9 K. However, this model was found to be thermally unstable due to the 
temperature-density dependence of the cooling rate by local bremsstrahlung emission. A related branch of hot optically thin 
disc solutions was found, when radiative cooling channels were extended towards higher temperature domains, including 
synchrotron radiation and comptonization of synchrotron and bremsstrahlung photons. These solutions are thermally and 
viscously stable for long-wavelength perturbations. The first models of this class were local thermal equilibrium models, i.e. the 
equation of energy balance is solved locally for any radius from a given underlying disc structure. Narayan & Yi (1994, 1995a, 
1995b) use a self-similar disc structure with sub-Keplerian angular momentum distribution (Begelman & Meier 1982; Spruit et 
al. 1987), whereas Abramowicz et al. ( 1995[ ) and Chen (199E) assume a slim disc structure with Keplerian angular momentum 
distribution. In spite of these kinematic and other minor differences in the particular cooling functions, both approaches yield 
the same local thermal equilibria ( Chen et al.1995] ). Advection becomes relevant if temperatures are very high, such that the 
relative half-thickness of the disc, h/r, is not negligibly small. This is the case in hot optically thin disc flows, where possible 
heating mechanisms primarily heat the ions to nearly virial temperature. Since Coulomb collisions become inefficient the 



e-mail: jpeitzaisw.uni-heidelberg.de 
t e-mail: appl@astro.u-strasbg.fr 



© 0000 RAS 



2 J. Peitz and S. Appl 



electrons gain only little of the thermal energy of the ions and remain at temperatures of T c ~ 10 9 — 10 10 K due to inefficient 
cooling. The local thermal equilibrium solutions are in general globally advection-dominated. Morphologically they resemble 
the ion-supported torus model of Rees et al. (1982). It turns out ( Abramowicz ct a!.199E| ) that the branch of stable optically 
thin equilibrium solutions is limited to a maximum accretion rate M max ~ a 2 , as long as the viscosity parameter a is smaller 
than some critical value a C rit, depending on the details of the model (Chen et al.1995). 

Recently, global solutions of hot optically thin accretion discs have been calculated. This requires to solve the full system 
of coupled equations of motion for the fluid in the disc, causing considerably higher numerical efforts due to the sonic point 
in the radial Euler equation. Viscous transonic solutions were calculated within the pseudo-Newtonian potential (Paczyhski 
& Wital980) for wedge-shaped flows by Chakrabarti (1996) and for disc-like flows by Chen et al. (1996) and Narayan et 
al. (1996). It turned out, that for most radii, the local thermal equilibrium solutions yield a good approximation for the 
temperature distribution in the disc. At the inner edge of the accretion disc, located at the horizon of the black hole, matter is 
forced to corotate with the black hole surface. This condition can not be implemented properly within the pseudo-Newtonian 
potential formalism. Therefore, disc properties in the very vicinity of the black hole are to be studied in full relativity. 
Furthermore both, stellar mass black hole candidates in galactic binary systems as well as the putative supermassive black 
holes in the centers of active galactic nuclei are expected to generally possess a considerable amount of angular momentum. In 
the first case the hole is spun up by accretion through the disc, whereas in the latter case the existence of jets gives theoretical 
support for a rapidly rotating central object (e.g. Blandford 1977; Rees et al. 1982). 

Accretion flows around rotating black holes have to be studied in Kerr geometry. A set of equations for stationary disc- 
like viscous transonic accretion flows around rotating black holes has been proposed by Lasota (1994; hereafter L94). Global 
transonic solutions of this system, corrected for some misprints, were given by Peitz (1994) for a polytropic equation of state. 
Solutions of a related set of equations, describing hot optically thin discs cooled by bremsstrahlung emission, have recently 
been calculated by Abramowicz et al. (1996, hereafter ACGL96) and discussed by Lasota & Abramowicz (1996). So far, most 
formulations for angular momentum transport have been using approximations for the shear, which are inconsistent with the 
kinematics in the supersonic inner part of the accretion discs. 

In this paper we derive the system of equations governing stationary viscous transonic disc flows in Kerr geometry 



(sections 



In particular, we give the consistent expression for the equation of angular momentum balance, important 



for the boundary condition at the horizon. We then discuss in section Q the properties of critical points in inviscid transonic 
disc-like flows, which hold as well for viscous transonic flows in the limit of small viscosity. In section ^ we describe a numerical 
method to solve the structure equations for transonic flows and present global solutions for a polytropic equation of state. 



2 BASIC EQUATIONS 
2.1 Metric 

The half-thickness h of the disc is assumed to satisfy the condition h/r 1. Then a power series expansion of the Kerr metric 
written in cylindrical coordinates {t,(j>,r, z}, to zeroth order in the small quantity z/r, yields a first approximation for the 
underlying disc metric. In geometrical units (G — c = 1), which are used throughout this paper, this metric g reads (NT73) 

g = - 1 —^dt 2 + ^-(d(f) - Ludt) 2 + ^-dr 2 + dz 2 , (1) 
A r 1 A 

where 

A = r 2 - 2Mr + a 2 , 

A = r 4 + r 2 a 2 + 2Mra 2 , 

v = -9t$/g<f>4 = 2Mar/A . 

Here M is the black hole mass and a — J/M its specific angular momentum. A vanishes at the event horizon, given by 
rh = M + (Al 2 — a 2 ) 1 / 2 . A subscript 'h' refers to quantities evaluated at rv 



2.2 Conservation laws 

Primary variables are the particle current n, the stress-energy tensor T and the entropy flux s. The motion of the fluid is 
governed by conservation of n and T, written as 

= V • n , (2) 
= h-(V-T), (3) 
= u-(V-T), (4) 
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where u is the four velocity of the fluid and h = u® u + g projects into its local rest frame (LRF). Equations (^), (^), (^) arc 
the continuity equation, the Navier-Stokes equations and the energy equation, respectively. 

The particle current is n = mi with the baryon number density n. The stress-energy tensor for a single component fluid 
with viscosity and heat flux is 

T = (p + p)u ®u+pg + t + q®u + u®q. (5) 

Here p — po + e is the total energy density, po = nms is the baryon rest mass density (ttib being the mean baryon rest 
mass), e is the internal energy density and p the isotropic pressure. All thermodynamic quantities refer to the LRF. Viscous 
contributions to T are contained in the viscous stress tensor t, given by 

t = -CO/I - 2ncr , (6) 

where a is the shear tensor, the scalar functions n and £ are the shear and bulk viscosity coefficients, respectively, and 
O = V • u is the expansion of the fluid world lines. The energy flux relative to the LRF is described by the spacelike vector 
q. The four-entropy flux is defined as 

s = nsu + q/T , (7) 

with s the specific entropy density and T the thermodynamic temperature. 

Stationarity and axial symmetry of g imply the existence of two commuting Killing vector fields, k = dt and m = d^,, 
which give rise to two conserved currents, 

P = T -k , V -P = , (8) 
J = Tm, V J = 0. (9) 

For ideal flows (t = q = 0), Eqs. (^|), (g) yield two constants of motion, namely energy E = —pu t and angular momentum 
L — (jmj,, where u t is the binding energy and p = (p + p)/po the specific relativistic enthalpy. The specific angular momentum 
A = L/E = —u^/ut is introduced for later use. 

2.3 Velocity field and reference frames 

We consider accretion flows with a four velocity field u — (u*, f2u*, u r , 0), where = u^/u* = d<f>/dt is the angular velocity 
and u r = dr/dr the radial velocity component, negative for accretion. In addition to the LRF attached to the fluid in motion 
with u, two other reference frames are frequently used in relativistic accretion disc theory. 

The locally non-rotating frame (LNRF) is characterized by u$ = and orbiting with angular velocity uj. Tensor indices 
with hat ' ' ' refer to LNRF. The physical 3-velocity with respect to LNRF has azimuthal and radial components 

l$) AO. if.) u T . . 

where Q = fi - u and U 2 = h rr = (u) 2 + A/r 2 . Here 

is the azimuthal Lorentz factor w.r.t. LNRF. Using 7 allows to write the normalization u ■ u — — 1 as 
t jA^ 2 U -yA^UQ 

u = A , UJ, = r2A ■ (12) 
In the following, we frequently write u = u r and l = u$. In terms of these variables, 7 and Cl are 

/ £^j2 \ !/ 2 r 2 Al 
1= \ l + MT>) ' Cl= A(AW + AP)^ ■ (13) 

The corotating frame (CRF) is orbiting with the angular velocity il of the fluid, being related to the LNRF by a Lorentz 
boost with 7 in the azimuthal direction. Tensor indices in CRF are labelled by bar ' ~ '. The radial component of the physical 
3-velocity w.r.t. CRF is = ^v {f) = u/U. 

The dragging of inertial frames forces any matter at the horizon to corotate with the surface of the black hole, Oh = or 
7h = 1, which follows from Eq. ([l3|). Consequently, v^ r \ and v^ r \ are both equal to the speed of light by virtue of the metric. 
This is in contrast to accretion onto non-collapsed objects, where corotation with its solid surface is imposed as a boundary 
condition. For later reference we further introduce the total Lorentz factor w.r.t. LNRF, 

7t = (l - (^) 2 - (v^) 2 ) ~ 1/2 = 7(1 - (^)T 1/2 , (14) 
which is singular at the horizon due to ?/ r 'h = 1. 
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3 EQUATIONS OF DISC STRUCTURE 

The flow is treated as a 1-D stationary problem, where the disc midplane coincides with the equatorial plane of the hole, mass 
increase of the hole and self-gravity of the disc are neglected. We use vertically integrated thermodynamic variables. Then 
the equations governing radial and vertical disc structure decouple. 



3.1 Continuity equation 

From the continuity equation (nu a )- a = it follows that the total flux of particles through a surface of constant r is a 
constant, the accretion rate M, given by 

M = -2-KruZJo , (15) 
where 

£o= Podz ~ 2hp (16) 

J-h 

is the surface rest mass density. For later use we introduce two other vertically integrated thermodynamic quantities, namely 
the total surface density, E ~ 2hp, and the vertically integrated pressure, P ~ 2hp. 

3.2 Radial momentum equation 

The radial Navier-Stokes equation K" a T a / — can be written as 

u'u + u 2 g+u 2 n = -(p+p)~ 1 (u 2 p' + s r + Q r ) , (17) 

with Q, Tt, K defined by 

g = 0nV5^), r = {a - Mr)/{rA) , (18) 
H = 7 2 AM/C/(r 2 A) 2 , (19) 

k = (n-ni){n-^)/(n+^). (20) 

The angular velocities 

n±_± Ml/2 (21) 

-± r 3/2 ±aM l/2 W 

correspond to prograde (+) and retrograde (— ) Keplerian orbits. Partial derivatives w.r.t. r are denoted by a prime. The 
quantity K measures the deviation of the angular velocity from the Keplerian value, with limiting cases /C(fi = nf) = and 
/C(H = 0) = 1. The quantities S r = h ra t a ^ and Q r — h ra (q a u /s + q^u )-^ describe contributions due to viscosity and heat 
flux. In the following we assume these forces to be small compared to pressure forces. Then the radial Navier-Stokes equation 
( p"7j ) reduces to the radial Euler equation for a perfect fluid. With the help of the continuity equation, and after integration 
over the disc height, the radial Euler equation can be be written as 

- = ^r> A/ = —y - H , T> = — - c ad , (22) 

u V r U 2 u 2 

where c a d is the adiabatic sound speed. Eq. (^) possesses a critical point for V = 0, identified with a sonic point in the CRF. 
Since causality requires c a d < 1, the conditions |i/ r 'hi = 1 and |f/ r '| < c a d at r ou t suffice to conclude, that black hole accretion 



is transonic (Shapiro & Teukolskyl983). Thus any global solution must satisfy 

Af{V = 0) = (23) 
as a regularity condition. 

3.3 Conservation of angular momentum 

The azimuthal Navier-Stokes equation can be expressed in the form of angular momentum conservation J.™ = 0, where the 
conserved current J r is given by 

J r = (p + p)u r U4, + + q U0 + q^u . (24) 

Neglecting the angular momentum associated with the heat flux and using the continuity equation yields 

= poru (put,) , r + (rt£), r . (25) 
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The vertically integrated angular momentum balance follows after integration over r, 

^(L - Lo) = 2hrtl , (26) 

where Lo is an integration constant. We did not neglect angular momentum associated with the internal energy. Therefore Eq. 
( f^l ) is the general equation of angular momentum balance (apart from a contribution due to photons) for any temperature, 
not restricted to flows with p ~ 1 as in NT73 or L94. 

3.4 Description of shear tensor and viscosity 

We assume the coefficient of bulk viscosity to vanish (£ = 0). The viscous stress tensor then reduces to t = —2r/cr. The shear 



tensor component <r£ can be written as (Eq. (A4)) 



2a\ = U 2 l (I + H\ + ^ - \ + 3G) + 2 [(Ma 2 - , 3 )fi - Ma] 2^ . (27) 
As demonstrated in App. [a|, expression (^) may be approximated by 

2cr; = U 2 l' + 2 [(Ma 2 - r 3 )fi - Ma] jA i/2 Ur-' L . (28) 

Based on the dimensional estimate for the viscosity due to the motion of turbulent fluid elements in the disc, 77 ~ po^turiAurb, 
we use a parametrization for the kinematic viscosity coefficient v = n/po, given by 

v = ac^U , (29) 

where a is the viscosity parameter and l tU rb ~ r 13 allows to mimic the standard assumption Z t urt> ~ h as well as other 



dependencies, such as Zturb ~ P/P' (Papaloizou & Stanleyl986). Viscosity is a transport process and operates on a finite 



time-scale, t„, limited by causality, no matter whether it is mediated by turbulent gas motion, magnetic fields, or photons. 
r„ becomes equal or larger than the accretion time-scale r acc as the accreted material approaches the horizon, and as a 
consequence the viscosity vanishes there, 

= . (30) 



A cut-off function f c has been introduced to guarantee such a behaviour. We choose (Narayanf992) 



/c = I (l-(« w /c) a ) ; \v^\<c v , (31) 
1 ; |t!< f >| >c v , 

parametrized by the maximum velocity of viscous information transport, c v . Since — 0, viscous stresses vanish at the 
horizon (no-torque condition; Abramowicz & Prasanna 1990) and the integration constant Lo in Eq. (^6|) is the angular 
momentum accreted by the black hole. Its value is obtained as part of the solution. 

In the limiting case of a Keplerian disc, i.e. for r > r ms and u r — > 0, vT r — ^ 0, i — > ik, — ^ ^k, Eqs. (p7|), (S) both reduce 
to the expression of NT73, 

2<t£ = 7k 3 J 4 3/2 A 1/ V- 5 ^ , (32) 

where 7k = 7(f2k) denotes the Keplerian Lorentz factor w.r.t. LNRF. For non-Keplerian accretion flows, L94 proposed to 
replace 7^ in Eq. ( |32| ) by the azimuthal Lorentz factor w.r.t. LNRF, 7, as in Eq. (jll|) and Q' k by Q' . Since 7h = 1, the resulting 
expression vanishes at the horizon, and thus leads to a no-torque condition as well. It does so by constraining the kinematics 
instead of using the physically motivated condition (|30|). In addition, this expression leads to a reversal of shear in the vicinity 
of the black hole, where 51' changes sign. Recently, ACGL96 replaced 7 in the L94 expression by the total Lorentz factor w.r.t. 
LNRF, 7 t, as in Eq. (|l4|). Different expressions for crj are compared in App. [a|. 



3.5 Conservation of energy 

a r 

L u;f3 



The energy conservation equation U a T^. s = can be combined with the continuity equation to 



p Q u r Ts, r = -o- a pt al3 - 1/362" - q a a a - q"* a , (33) 

where a = V„u is the four-acceleration of the fluid. The compressive heating rate —1/30^ vanishes due to the assumption 
C = 0. The special relativistic term a ■ q describes the inertia of the heat flux. Neglecting this term, the vertically integrated 
energy equation is given by 

Qadv = Q+ - Q- , (34) 
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where Q a d v = uSoTs tT describes the advection of entropy by the radial motion of the fluid, Q_| 
viscous heat generation and Q- = —2hq z — 2hr~ 1 (rq r )^ r is the local cooling rate. 



2ha a pt al3 is the rate of 



3.6 Equation of state 

As a first approximation we assume a polytropic relation between the vertically integrated pressure P and the surface rest 



mass density So (Popham & Narayanl991) 



2 
f -'ad 



P = KEq , 

where K(s) is a constant measuring the entropy of the flow and the polytropic index for a two-dimensional flow, N = l/(-T 
is constant. The adiabatic sound speed may then be written as c ad = (dP/dS) a d, leading to 

t dP \ s _ rp/s 

\dS ) ad S + P l+NPP/So' 
Using the continuity equation then yields 

The 'modified accretion rate' 
M =r N K N (M/2n) 



(35) 
-1), 



(36) 



(37) 



consists of a combination of constants M, K and r and thus is a constant. For fixed M and P, A4 measures the entropy in 
the flow. For fixed K and P, M is a measure for the accretion rate. For P > 1, the specific relativistic enthalpy in terms of 



c ad is /* = (1 - iVc a 



The polytropic equation of state directly couples the temperature to the kinematics of the flow, and 



makes the energy equation obsolete. Instead of two independent parameters M and K, one is left with a single parameter M. 



3.7 Vertical disc structure 

The fluid is assumed to be in vertical hydrostatic equilibrium, i.e. in the LRF the vertical pressure gradient p tZ balances the 
vertical tidal acceleration towards the midplain of the disc. A proper treatment of this force balan ce requires an expansion 
of the z-component of the Navier-Stokes equation to O(zfr) 1 . As shown by Riffert & Herold ( 1995 ), in a Keplerian disc the 
hydrostatic equilibrium, obtained from the vertical Euler equation in this way, differs considerably from the expression given 
by NT73. For general disc flows with u = (it*, , u r , 0) the vertical Euler equation ^] h % T!jf = 0, expanded to O(zfr) 1 , may 
be written as (App. H) 



1 



dp 
p + p dz 

PL 



2 7 fAA 1 



y = i- 



[y(l - aQ) 2 - 4af7(l 
2M 3a 2 



a,Q) + 2r 2 Q, 2 ' 



A^U 2 



Z = 



2M 
r 



(38) 
(39) 
(40) 



It can be shown that TZ remains finite at rh and therefore (B8J)is regular everywhere. In the limiting case of a (cold) Keplerian 
disc flow, e.g. for c a( j —> 0, u — > 0, Q. — > f2k, p + p — > po, Eq. (|38|) reduces to the formula given by Riffert & Herold (199E). To 
obtain an explicit expression for the half-thickness h from Eq. (^), we substitute dp/dz ~ p/h and assume 

S P S t dP\ 



V 



p + p S + P S 
This leads to 
h 2 



S + P \dS0J; 



_ 2 

— c ad 



M1l + cl d r 2 Z(Af/V + l/r) 



(41) 



(42) 



If the radial pressure gradient in the vertical Euler equation is negligible, (^2|) reduces to h 2 /r 2 — (r/A/)c ad ^.~ 1 . 

The incomplete NT73 expression for the vertical hydrostatic equilibrium in Keplerian flows, given for comparison with 
alternative expressions for the disc height, reads 

(r 2 +a 2 ) 2 + 2Aa 2 ' 



po dz r r 2 



ft = 7 k 



(r 2 + a 2 ) 2 



Aa 2 



(43) 



t The derivation of the thickness of transonic accretion discs form the verti cal Euler equation was proposed by J. P. Lasota, who also 
drew our attention towards a different approach by Abramowicz et al. (1996). 
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where 7k is the Keplerian Lorentz factor w.r.t. LNRF. For non-Keplerian accretion discs, L94 proposed to replace 7^ in Eq. 
( p3[ ) by the azimuthal Lorentz factor (p"l|). Recently, ACGL96 used the total Lorentz factor ( |l4| ) instead of 7^. Since 7t is 
singular at the horizon, the disc half-thickness h, as derived with this approximation, vanishes at r^, which leads to infinite 
densities and therefore to difficulties in the energy equation. The relative disk height (^2|), as well as alternative expressions, 
are shown in Fig. 12. Equation (^2|) yields the complete expression as long as the vertical velocity component u z can be 
neglected. 



4 SOLUTION TOPOLOGIES 



The regularity condition at the sonic point allows to identify different solution topologies for accretion flows (Abramowicz & 



Zurekl981). This analysis, performed for a pseudo-Newtonian potential, as well as related studies for wedge-shaped flows in 



Kerr geometry (Chakrabartil990) are generalized in this section. A subscript 's' refers to 'critical quantities', i.e. quantities 
at r s . 

As shown in App. |c|, the regularity condition Af(V = 0) = selects a 2-D hypersurface T of the parameter space spanned 
by three critical variables, e.g. by (r s , c s , Q s ) or by (r s , E s , A s ). Projections of T onto the E s -X s plane are plotted in Fig. 1 
for a non-rotating and for a rotating (a — 0.95) black hole. For E s , A s within ABC, three solutions for r s exist. These are 
referred to as inner critical point, r s ,i n , intermediate critical point, r B ,int, and outer critical point, r s , out , according to their 
relative location. For E s , A s beyond the lines BC and BA, only one solution for r a exists, corresponding to an inner critical 
point for E a , A s within / and to an outer critical point for E s , A s within O. Inner critical points are located in the vicinity of 
the marginally stable orbit r ms (r ms = 6r g for a = and r ms ~ 1.94 r g for a — 0.95, r g = M being the gravitational radius), 
whereas outer critical points are located farther away from the black hole (App. |c|). The type of critical points depends on 
the critical slope (dX/dr) s , and separates unphysical from possible physical solutions. 

In ideal flows, with A = A s and E = E s being conserved, r s ,j nt is an center-type critical point. Any other critical points 
are saddle-type, corresponding to possible physical solutions. Ideal solutions are determined by two out of three constants of 
motion E, A, A4. For given E and A, the corresponding accretion rate M increases as r s decreases. For any A within ABC, 
there exists an energy E m such that the corresponding A4 for a solution passing r s ,j n is equal to A4 for a solution passing 



r 8 ,out) -M m (i5 m , A, r S] i n ) = Mm(E m , A, r Sj0 ut) (Chakrabartil99C ). Thus the corresponding line BD indicates, where stationary 



solutions actually 'jump' between different topologies as E, A are varied continously. In the following, solutions with E s > E m 
are referred to as type I , and solutions with E s < E m as type II . Solutions with E s = E m pass through both r Sj i n and r BjOU t 
and are therefore global accretion solutions only if a shock occurs within r s ,in and r SjOU t. The dashed line in Fig. 1 indicates 
the maximum A which a global ideal transonic flow with a given E may have in order to pass the centrifugal barrier. 

Viscous solutions with non-conserved E and A are represented by curves E(X). The location of r s results from the position 
of i?a and A s on that curve, which is generally unknown in advance. The solution topology at r s , depending on (d\/dr) s , is now 
determined by the equation of angular momentum balance, introducing additional parameters, e.g. a and Lo- For moderate 
(dX/dr) s , the topological classification in terms of local critical quantities E s and A s yields essentially the same result as for 
ideal flows. The E(X) curves of such flows may then be superposed on Fig. 1. If E(X) never crosses the line BD, the solution 
topology and the approximate location of r s are known. This situation applies e.g. for quasi-ideal flows with globally small 
gradients dE/dr and dX/dr, where E(X) is bound in narrow regions around E s and A s . If E(X) crosses the line BD, the 
unknown position of E s and A s becomes relevant. In special cases E s and A s are restricted, e.g. if the flow is assumed ideal 
within r s . 

Figure 1 shows that flows around rotating black holes possess generally higher energies and thus higher temperatures 
than flows around non-rotating black holes with corresponding topology. 



(44) 



5 GLOBAL SOLUTIONS 

5.1 Model equations and numerical method 

The system of two coupled ODEs solved in this section consists of Eq. (J^^ and Eq. (^6|) combined with Eq. (Pq), 

du uN 
= "iT ' 

Tr - W^^^-M^]^. ,45, 

The adiabatic sound speed c a d is calculated according to Eq. (^tj), assuming iV = 3/2 (T — 5/3). The kinematic viscosity 
coefficient v is parametrized by Eq. (^), using c v = 1 and (3 = 1. The inner edge of the disc is located at n n = r^, the outer 
edge at r out = 100r g . The radial Euler equation has a sonic point at radius r s , which is calculated as part of the solution by 
the condition T> = 0. This provides, together with the regularity condition N(D = 0) = 0, the boundary condition for Eq. 
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(p^). At rin the condition = requires Lh = Lo, which is used as a boundary condition. The four boundary conditions for 
the four unknowns it, I, r s , Lo are therefore given by V = 0, M — 0, L h — Lo and e.g. A(r ou t) = A ou t- 

Eqs. (^), ([IB]) are mapped onto the unit interval in x, separately for the sub -and supersonic domain. The resulting four 
equations are then solved with x as an independent variable. The sonic point is mapped onto x = 0, where M = 0, T> = are 
taken as boundary conditions both for the supersonic and subsonic equations. The horizon and the outer radius of the disc 
are both mapped onto x = 1, where L — Lo is used as boundary condition for the supersonic equations. For a given Lo, the 
solution is then uniquely determined, including A ou t- An iteration scheme in Lo is set up to find global solutions for a given 
Aout- The boundary value problem is integrated by a collocation method (Bader & Kunkell989), using a fully adaptive grid 
with variable number of subintervals. This allows to obtain the solution to a pre-specified accuracy. A parameter continuation 
is used to study the variation in a, a and Ml, as well as for the iteration in Lo- 

We use dimensionless quantities by normalizing distances to the gravitational radius r g and velocities to the speed of 
light (geometrical units). For a polytropic equation of state the resulting system is independent of mass M. We present ideal 
and viscous solutions for any of the topologies discussed in Sect. ^[ The relative half-thickness h/r of the disc is calculated 
from Eq. 042). No retrograde disc flows with A < are considered. 



5.2 Ideal solutions 

To give examples for all solution topologies, we plot model sequences of increasing A for E — 1.015 in the case of a non-rotating 
(a = 0) black hole (see Fig. 2) and for E — 1.06 in the case of a rotating (a — 0.95) black hole (see Fig. 3). The various 
models correspond to crosses x in Fig. 1. For lower A, within ABD or O, solutions pass an outer critical point. For A within 
O, u grows with decreasing r. For A within ABD, the existence of the inner and inter sonic point becomes 'visible' in a 
deformation of u at the corresponding locations. As A approaches the line BD from within ABD, u passes the outer sonic 
point followed by a local maximum and a local minimum, both within the supersonic region. This behavior occurs because 
M = is encountered in total three times, first together with T> = at the sonic point, and thereafter two more times for 
finite T>. As A crosses the line BD, stationary solutions 'jump' from type II to type I topology, passing the inner critical point 
(compare e.g. Chakrabarti 1990; Kafatos & Yang 1994). As A approaches the line BD from within BCD, u first passes a local 
maximum followed by a local minimum, both within the subsonic region, before passing the inner sonic point. The majority 
of type I solutions with A within / have again monotonic u-profiles. 

The 'jump' of stationary solutions at the line BD is not visible in the corresponding S7-profiles. In the vicinity of the 
horizon, the ^-distribution shows a kinematic boundary layer structure, dropping with a steep gradient to Qh = 0. In that 
kinematic boundary layer, the maximum of Q and the steepness of dQ/dr depend on the angular momentum of the flow and 
the rotation of the hole. This is illustrated in Fig. 4, where £1 is plotted for a sequence of type I models with increasing a and 
fixed E = 1.2, A = 1.0. 



5.3 Viscous solutions 

For given boundary conditions the global solution is uniquely determined and possesses either type I or type II topology. 
Type I solutions are generally found for moderate a and higher M. In Fig. 5 we plot type I solutions around a non-rotating 
black hole, which differ only in a, with otherwise fixed a, Mi and Lo- They show a globally sub-Keplerian angular momentum 
distribution. With increasing a the angular momentum transport becomes more efficient, such that practically all reasonable 
boundary conditions A ou t, from A ou t <; A s to A ou t ^ Ak(r out ) are met. For a fixed, A ou t increases with Lo- This is demonstrated 
in Fig. 6 for a sequence of type I solutions around a rotating (a — 0.95) black hole, which are for the same a, a and M, 
but differ in Lo. There exist solutions which approach a Keplerian distribution Ak globally or even exceed Ak locally. These 
require high A in the vicinity of r ms , obtained for lower a and possess higher Lo. In Fig. 7 we plot two examples of such 
trans-Keplerian solutions for a = 0. Solutions of that kind possess extremely high temperature at r out , such that their physical 
relevance appears doubtful. For very low viscosity, e.g. for a = 10 -4 — 10 -3 , accretion profiles are quasi- ideal. Under these 
conditions, type I solutions with E(X) within BCD can be calculated, possibly passing local extrema in their it-profiles. In 
particular, the quasi-constant A may exceed Ak(r ms ) in the vicinity of r ms . In Fig. 8 we plot examples of type I quasi-ideal 
solutions for a — 0, with globally sub-Keplerian as well as locally super-Keplerian A-distribution. For higher A or for energies 
below unity, r s ,out is shifted outwards, in extreme cases beyond r out (Fig. 14). Then even solutions with inverted angular 



momentum transport in their sub-Keplerian outer parts exist (Fig. 8), as reported by Abramowicz et al. (1988). However, in 
order to guarantee that these are global solutions, the integration must be extended to larger r out . 

Except from quantitative differences, type I solutions show the same qualitative behavior for any a. The Cad-distribution 
and thus the temperature distribution possess a maximum outside r s . For given a, M and Lo, higher a yield steeper A-profiles 
and higher temperatures. Therefore, quasi-ideal flows represent the limit of lowest temperature discs for a set of a, M and 
Lo. Since temperature remains high at r out , no (cold) standard thin disc configuration is met. This reproduces a well-known 
property of hot accretion flows, and is not an artefact of a polytropic equation of state. A common feature of type I solutions 
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is a very flat A-distribution inside r s . This is due to low v resulting from the decrease of c a d within r s . The energy in the 
disc increases moderately towards r out . The decrease of c a( j towards is imprinted on the E- profiles and the //-profiles (not 
plotted), which show steeper gradients than A in the vicinity of r^. Thus, in general, even for quasi- flat A-distribution inside 
r s , the angular momentum Lo deposited onto the black hole in a type I flow is significantly lower than L s . The relative disc 
height respects h/r & 1 at r s , but in models with effective angular momentum transport exceeds unity beyond ~ 20r g due to 
very high temperaturs. The ^-distributions show a kinematic boundary layer structure similar to ideal flows as a consequence 
of flat A-distributions within r s , such that disc models with v = in the supersonic domain (Popham & Narayanl992), 
essentially yield similar profiles. 

Viscous solutions passing the outer sonic point exist for higher a and lower M . Since the angular momentum at r SjOU t 
is strongly sub-Keplerian, type II solutions possess smaller Lo. In Fig. 9 we plot a sequence of type II solutions around a 
non-rotating black hole, for fixed A4 and with different Lq. Since r s 2> r^, global solutions with any reasonable boundary 
condition A out do not exist for each a. The two solutions with higher A in Fig. 9 are for the same a = 0.06, choosen such 
that the corresponding maximum possible A ou t is Keplerian. Lower a only allow for A out < Ak(r ou t) (for otherwise fixed 
parameters). The two solutions with lower A in Fig. 9 are for a = 0.05. The global A = distribution with Lo = does 
no more depend on a, describing the limit of 'spherical disc accretion' around a non-rotating hole. This limit is reproduced 
for any E and A within Fig. 1. Alternatively, as for type I solutions, different A out can be found for fixed a, M and Lo by 
variation of a. This requires A4 high enough in order to achieve r s < r ou t- A sequence of such type II solutions around a 
rotating (a = 0.95) black hole is shown in Fig. 10. As for type I flows, quasi-ideal solutions with E(X) within ABD exist, 
again similar to corresponding ideal solutions and not shown here. 

Solutions around rotating and non-rotating black holes differ only quanitatively. Since the A-distribution is strongly sub- 
Keplerian at r s , accretion flows with smaller A out result in general and no type II solution with a global Ak-distribution does 
exist. Nevertheless, A can be steeper within r s than for type I solutions due to higher a. Type II solutions have generally 
lower temperature than type I , with a maximum at the horizon. However, as before, no (cold) standard disc configuration 
is matched at r ou t- With decreasing M, r s moves outward and c a d decreases due to lower A s and thus less efficient global 
angular momentum transport. Therefore the relative disc height h/r exceeds unity only in the very outer disc regions. The 
kinematic boundary layer is less prominent due to generally lower A. We do not find solutions with negative accreted angular 
momentum. 



6 SUMMARY AND DISCUSSION 

We derived the equations of motion for stationary viscous disc-like accretion flows around rotating black holes. In particular, 
we give a complete expression for the shear tensor that governs the transport of angular momentum in relativistic transonic 
disc flows. In order not to violate causality, the kinematic viscosity has to vanish at the horizon. This is achieved by an 
appropriate parametrization. Consequently, no torque is exerted on the horizon. It is shown that the expression widely used 
for the shear, calculated by assuming circular trajectories, is incompatible with the kinematics in the supersonic inner part of 
the accretion discs. We wish to stress that it is important to use the correct boundary conditions at the horizon, even if it turns 
out that they do not strongly affect the global properties of solutions. The range of parameters, for which global solutions 
exist, can be affected by the choice of boundary conditions. Furthermore, we calculate an expression for the half-thickness of 
the disc, obtained by integration of the vertical Euler equation over disc height. 

The properties of the sonic point are exploited to classify different solution topologies of ideal disc-like accretion flows 
around rotating black holes. Depending on angular momentum and energy, global solutions pass through the inner or outer 
critical point, referred to as type I and type II , respectively. We use a polytropic equation of state for vertically integrated 
pressure and density to calculate global viscous solutions. The coupled disc equations are solved between the horizon of the 
black hole and 100r g . A numerical method is presented that allows to integrate the structure equations of transonic accretion 
flows. In cases of moderate angular momentum transport, the topological classification for ideal flows also holds for viscous 
flows. We give examples of each topology for non-rotating and rapidly rotating black holes (a = 0.95). We present solutions, 
where the angular momentum at the outer edge of the disc varies between zero and the Keplerian value, for a wide range 
of parameters. The angular momentum Lo accreted by the black hole is calculated as part of the solution. For small values 
of a and higher M, solutions are of type I , with the sonic point located in the vicinity of the marginally stable orbit. In 
these cases, the temperature distribution in the disc has a maximum outside the sonic point. There exist solutions which 
are sub-Keplerian everywhere (Figs. 5, 6), and accretion flows, which are super-Keplerian in some region inside and outside 
the marginally stable orbit (Fig. 7, 8). In this case, the angular momentum accreted by the black hole generally exceeds the 
Keplerian value at the marginally stable orbit. For these partially super-Keplerian flows the sonic point is located between the 
photon orbit and the marginally stable orbit. For a given angular momentum at the outer edge of the disc, trans-Keplerian 
flows possess a higher temperature than sub-Keplerian flows. For higher values of a and smaller M, solutions are of type II 
, with the sonic point located further away from the black hole. They generally possess a sub-Keplerian angular momentum 

© 0000 RAS, MNRAS 000, 000-000 



10 J. Peitz and S. Appl 



distribution for any radius. These solutions are colder than type I and have their maximum temperature at the horizon. 
In the outermost part of type I solutions, the relative disc height h/r exceeds unity, or the temperature reaches the virial 
temperature. Our solutions however remain consistent for the inner 20r g . 

We obtain solutions around rotating black holes for any Kerr parameter a. Discs around rotating black holes are hotter 
and deposit less angular momentum onto the central object than discs around non-rotating black holes. For the latter the 
temperature maximum is broader and at larger radii than in the case of a rapidly rotating black hole. Temperatures of our 
solutions lie in the range of 10 11 - 10 12 K. They remain high throughout the whole disc and our solutions do not match 
a (cold) standard thin disc configuration at the outer edge. These results suggest to compare our polytropic model with 
other hot accretion disc models, e.g. with optically thin advective disc models. Global solutions of this type, described by a 
related system of equations, have been calculated recently (Abramowicz et al.1996). It turns out that the general features 
of both models argee. A polytropic equation of state therefore provides a reasonable description for many properties of hot 
optically thin advective disc flows. Since we do not find any solution with considerably lower temperature than ~ 10 11 K, 
a polytropic equation of state does not seem to reproduce the class of cold radiation-pressure supported optically thick slim 
discs ( Abramowicz et al.l988| ). 

We are grateful to M. Camenzind (who initiated this work) and to R. Khanna for invaluable discussions. Stimulating 
comments by S. Chakrabarti, J. Ferreira, J. P. Lasota and L. Rezzolla are acknowledged. This work has been partially 
supported by Deutsche Forschungsgemeinschaft (SFB 328) and the French Ministery of Foreign Affairs. 



APPENDIX A: SHEAR TENSOR COMPONENT <xf 



To calculate the shear tensor <j it is convenient to start from the decomposition of the gradient of the four- velocity, Vit, into 
irreducible tensorial parts, 



«a;/3 = &af3 + W a p + -@h a /3 — a a U/} , 



(Al) 



where w is the antisymmetric rotation 2-form, cr is the tracefree symmetric shear tensor and a is the four-acceleration of the 
fluid. The expansion of the fluid world lines is given by Q = u" a — vT r + u T jr. Symmetrization (expressed by brackets) of Eq. 
(IaI) leads to 



O a p = U( Q;(3 ) + a (Q U (3 ) — -@h a /3 . 

The relevant component for the angular momentum equation is then 

9 

r» r T . rr . r , r r 

la^ = u.^ + g u^-r + a + a^u — u <t> , 

which, for u = (u l ,flu ,u r ,0), can be written as 



I' 



u 2 l ( u' 



2a% = U 2 l [j+n + ^ NL_f+30 +2 [(Ma 2 -r 3 )fi- Ma] 



(A2) 



(A3) 



(A4) 



Instead of using this full expression (with u' /u = N /T>), crj may be approximated by a simplification of the velocity field u 
involved in the calculation of a^. The simplest approximation of a Keplerian disc flow with u r = leads to Eq. (^). While 
consistent for nearly circular trajectories given in Keplerian discs (NT73), Eq. ( |32j ) breaks down in the inner part of transonic 
discs, where radial and azimuthal velocities become comparable. 



A more appropriate simplification of Eq. (A4) consists of dropping terms which contain u' /u. We therefore neglect the 
expansion (O = 0) and the radial four-acceleration (a r — 0). This approximation yields 



2a% = U 2 l' + 2 [(Mo 2 



r 3 )n 



■ Ma] -/A 1/2 Ur' i 



(A5) 



We found that viscous solutions calculated using Eq. (A5) are consistent with the assumptions made above for all parameters 



and for almost any radius. In Fig. 11 the exact expression (A4) was used to calculate a typical type II solution. This solution 



was used to calculate the approximation (A5), as well as other expressions proposed for the shear. The exact expression (A4) 



is plotted as a solid line and re mai ns finite at the horizon. The approximation (A5) is plotted as a short dashed line and 
yields a good approximation to (A4) for all radii. Eq. (^) with 7 and fl' instead of 7k and fi k (long dashed line; L94) differs 
considerably, especially in the vicinity of the kinematic boundary layer. Eq. (ji^) with jt and fl' instead of 7k and fl' k (short 
dash-long dashed line; ACGL96) remains a good approximation for larger radii, but tends to extremely high values of inversed 
shear in the inner parts of the flow, due to the singular behaviour of 7 t at rh. For larger radii all expressions approach the 
Newtonian limit. 
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APPENDIX B: VERTICAL EULER EQUATION 



In this section the vertical Euler equation h%T^f = is derived for a velocity field u = («*, v?, u r , 0), to an accuracy of first 
order in (z/r). This requires an expansion of the Kerr metric to 0(z/r) 2 (Riffert & Heroldl995), denoted by G, and written 

as 

z z 2 ~ 
G = g + - g+^r g , 



where g is the zeroth order (Eq. 
Alu 2 Ar 2 



NT73), with non- vanishing components 



9tt = 



_ _ _ - _ (i _ 2M \ 



gtj, - gcj>t - --5- 



2Ma 



A 

= 3- 



A 



= 1 



g is the first order correction, with non-vanishing components 



Qrz — Qzr — Qrr 



and g is the second order correction, with non-vanishing components 
2a 2 " 



g tt = 

9t<j> = 

g^if, 

lirr 



r 



3 a 2 

9<t,t = ~9H I 9 + 72 



= -a 1 + 



2 

5M 2Ma 2 



r .3 



3M 4M 
r 



2 



4l3- 



6M 2a 2 
+ — 



2M 2Ma 2 a 



Neglecting contributions of viscosity and heat in T, the leading terms of h z a T a ^ — are 
G az p :a + (p + p)u a u z a = . 

The semicolon denotes covariant derivatives w.r.t. G. This yields 

G rz p, r + p, z = -(p + p) [T z t {u) 2 + r^uV + r^(ir ) 2 + rj; r (u r ) 2 ] , 

where we already used G zz = 1 + 0{z/r) 2 and rf r = F^ r = 0. The remaining expressions in Eq. ( |B6| ) are 

2r^,3 = — G a /3,z — G rZ G a f3,r , O, /? 6 {t, 4>} , 

G rz = -(z/r)£ + e>( 2 /r) 2 , 
which leads to 

r« = ~y + o(*/r) a , 



r? 



r 2 



a(2 + y)+C>(z/r) 2 , 
[a 2 (4 + ^) +2r 2 l +0(z/r) 2 



(Bl) 



(B2) 



(B3) 



(B4) 

(B5) 

(B6) 

(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 



where y, Z are given in Eq. (^). Equation ( |B(j ) is further reduced by = fiti* with u* as in Eq. (jl^) and the radial pressure 
gradient is, as in the case of the radial Euler equation (E2J) expressed in terms of c a d, 



1 dp 



2 

~ c ad 



— + - 

u r 



p + p dr 

This then leads to the final expression (M). The limiting value of TZ at the horizon, 1Z^, 



(l h - au h Zl /2 ) 2 - ih(l - Zh) 



4M 2 
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remains finite. 

In Fig. 12 we plot three different expressions proposed for the relative half-thickness h/r, calculated for a typical type II 
solution using the exact expression for the shear (same solution as in Fig. 11). The solid line corresponds to Eq. (fl2|), which is 
to be compared with the corresponding expressions derived from Eq. (ji^). These are given by (h/r) 2 = (r /M)c 2 d /lZ, with 1Z 
calculated using either 7 instead of 7k (short dashed line; L94) or 7t instead of 7^ (long dashed line; ACGL96) . All expressions 
yield similar flow morphology in the outer subsonic parts of the flow. The disc height (^) is finite at r n , in agreement with 
(B14). Using 7 in (0) overestimates h within ~ 10r g , whereas j t yields an underestimate, with vanishing h at the horizon. 



APPENDIX C: EXISTENCE AND LOCATION OF SONIC POINTS 

In the Kerr metric circular equatorial geodesic orbits exist outside of the photon orbit ri = 2M(1 + cos(2/3cos _1 (±a/M)). 
At r p h the Keplerian angular velocities O.^ reach the limiting angular velocities fl± for timelike trajectories (dashed lines in 
Fig. 13, bottom). Circular geodesic orbits are unstable within the marginally stable orbit r ms > r pn . In the following we use 
abreviations A ph = A k (r- pn ) and A ms = A k (r ms ). 

For given r s , c s the regularity condition N(D — 0) = yields a quadratic in fl s , 

= a + aifi s + a 2 fi? , (CI) 

a = c 2 r 2 (M - r s )(4AfV - r 2 A) - MA A , 

01 = 2MaA[A-2r s c 2 B (M -r s )] , 

a 2 = A[A(rf - Ma 2 ) + Ac 2 (M - r B )] . 



In the limiting case c s = 0, solutions of Eq. (CI) reduce to Slf 1 = fij, where fl+ > to corresponds to prograde and Q B < to 



to retrograde critical points. In the general case < c 2 < 1/3 one finds fij > fl+ for r£ > r^ h and fi k < fl s for r s > r ph 



Therefore transonic disc flows have sub-Keplerian f2-and A-distributions in the vicinity of the sonic point (Abramowicz & 



Zurekl981). A minimum critical radius r s _(c s ) exists, which is located between r s _ = r* h for hot flows with c 2 — 1/3 and 

0. The corresponding A represents an upper limit A top allowed for an ideal transonic 



polytropic flow with a given enthalpy (Kafatos & Yangl994). On the other hand, a maximum critical radius r s +(c s ) exists, 
which is smallest for hot flows and tends to infinity as c s vanishes. In Fig. 13 we plot fl*(r s ) and A*(r a ) for different values of 
c 2 as a parameter, both for a non-rotating and a rotating (a = 0.95) black hole. A limiting critical sound speed c*(a) separates 
the topology of the X s (r s ) curves in two classes. If c s > c*, A s (r s ) is monotonic, thus any A s < A p h corresponds to one solution 
for r s and A s > A p h yields no solution at all. On the other hand, if < c s < c*, A s (r s ) has a minimum A m i n and a maximum 
Amax, which fall together for c s = c*. Obviously, A m i n < A ms and A max > A ms . Consider first the case A max < A p h. Then for 
A s < A m i n again exactly one r s exists, for A m i n < A s < A max there exist three solutions for r s , and for A s > A max again one r s 
exists as long as A a < A p h, and no solution otherwise. The number of critical points in each of these three domains is reduced 
by one if A max > A p h- These critical points can be identified by the corresponding slope d\ s /dr a and their relative location. 
An outer critical point r StOUt exists under multiple critical point conditions and has d\ a /dr s < 0. An intermediate critical 
point r S] int has d\ s /dr s > 0. An inner critical point r Si i n has dX s /dr s < and is located next to the horizon. The simple 



relation (C2) among r s , c s , is locally relevant only. This changes for ideal flows if N(D = 0) = is solved for conserved 
quantities, since then the solution topology is known from the boundary conditions. Out of numerous combinations we plot 
^^(^a) parametrized by E s in Fig. 14 and E B (r a ) parametrized by A s in Fig. 15, both for o = and for a = 0.95. Again, the 
same classes of critical points are recovered from both diagrams. For E s < E* there exist extreme values A m i n and A max , such 
that for A m i n < A s < A max multiple critical points occur. Formulated vice versa, for A s > Ag there exist extreme values -E m i n 
and £ max , such that for i? m i n < E s < -E max multiple critical points occur. This is summerized in Fig. 1, where -E m i n (A s ) and 
-B max (A s ) are plotted. 

Consider now ideal flows with A and E in the region of three critical points, and ask through which of these a particular 
transonic solution passes. The answer depends partly on the nature of critical points, which results from the critical slope 
(du/dr) s , calculated by de l'Hopital's rule from critical quantities. It turns out, that r s ,i nt corresponds to an unphysical 
center-type critical point (spiral-type in viscous flows) , whereas the inner and outer critical points are saddle-type, allowing 
for physical solutions. Thus, as A increases continuously, the stationary solution 'jumps' from passing r s . out towards passing r Sj j n 
somewhere between A m i n and A max . Where this transition occurs in the _B-A-plane, depends on the modified accretion rate M 
corresponding to a particular r s . M is generally different for solutions passing r s ,i n or r Bj0u t, M(E, X, r s ,in) 7^ M(E, X, r s , out ). 
However, there exists an energy E m , such that A4 m (-B m , A, r s ,in) = Mm(S m , A, r SiOUt ). The line BD in Fig. 1 represents the 
locus of these E m . 
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APPENDIX D: FIGURES 

Figure 1: Division of the parameter space spanned by E s , A s according to different types of ideal flow, for black holes with 
a = (top diagram) and a — 0.95. Curve ABC is the boundary of the parameter space which seperates the region (between 
AB and CD) with multiple critical points from the region with a single critical point (rest of diagram). BD is the locus of 
flows for which both, E E and M, at the inner and at the outer critical points are the same. The dashed line indicates the 
maximum A which an ideal flow with a given E may have in order to pass the centrifugal barrier. The locations of ideal models 
presented in section H are indicated by x . 



Figure 2: Sequence of ideal disc solutions around a non-rotating black hole. Solutions are for a = and E = 1.015, with 
A increasing (from bottom to top in the lower diagram) as A = 3.0, 3.25, 3.33, 3.335, 3.5, 3.8. The sonic point is marked by 
• . The location of these models in the L s -A s -plane is marked by x in Fig. 1. 



Figure 3: Sequence of ideal disc solutions around a rotating black hole. Solutions are for a — 0.95 and E = 1.06, with A 
increasing (from bottom to top in the lower diagram) as A = 2.0, 2.05, 2.08, 2.09, 2.15, 2.3. The sonic point is marked by •. 
The location of these models in the _B a -A a -plane is marked by x in Fig. 1. 



Figure 4: Sequence of type I ideal disc solutions around rotating black holes. Solutions are for E = 1.2 and A = 2.0. The 
Kerr parameter a increases (from bottom to top) as a — 0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.95. The sonic point is marked by •. 



Figure 5: Sequence of type I viscous disc solutions around a non-rotating black hole. Solutions are for a — 0, M — 0.032 
and Lo = 3.5, with a decreasing (from top to bottom) as a = 0.045, 0.04, 0.03, 0.01. Dashed curves represent Keplerian 
distributions. The sonic point is marked by •. 



Figure 6: Sequence of type I viscous disc solutions around a rotating black hole. Solutions are for a = 0.95, Mi = 0.027 
and a = 0.022, with Lo decreasing (from top to bottom in the A-diagram) as Lo = 2.43, 2.40, 2.30. Dashed curves represent 
Keplerian distributions. The sonic point is marked by •. 



Figure 7: Sequence of type I viscous disc solutions around a non-rotating black hole. Solutions are for a — 0, M — 0.032 
and Lo = 5.0, with a decreasing (from top to bottom in the A-diagram) as a — 0.008, 0.006. Dashed curves represent Keplerian 
distributions. The sonic point is marked by •. 



Figure 8: Sequence of type I viscous disc solutions around a non-rotating black hole. Solutions are for a — 0, a — 0.001 
and Lo = 3.8, with M decreasing (from botom to top in the A-diagram) as M — 0.0112, 0.0028, 0.0010, 0.0004. Dashed curves 
represent Keplerian distributions. The sonic point is marked by •. 



Figure 9: Sequence of type II viscous disc solutions around a non-rotating black hole. Solutions are for a — and 
M = 0.011. The upper two solutions have a = 0.6 and (from top to bottom in the A-diagram) Lo = 1.3, 1.1. The lower 
two solutions have a = 0.5 and (from top to bottom in the A-diagram) Lo = 0.9, 0.0. Dashed curves represent Keplerian 
distributions. The sonic point is marked by •. 



Figure 10: Sequence of type II viscous disc solutions around a rotating black hole. Solutions are for a = 0.95, M = 0.011 
and Lo = 1.0, with a decreasing (from top to bottom in the A-diagram) as a — 0.8, 0.7, 0.3. Dashed curves represent Keplerian 
distributions. The sonic point is marked by •. 
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Figure 11: Different expressions for the shear tensor component <r£, plotted for the supersonic part of a typical type 
II solution (a — 0, a — 0.5, M = 0.017, Lo — 0.9). The solid line describes the full expression (A4), used to calculate 
the underlying disc solution. The short dashed line corresponds to approximation (A5). Other curves represent expressions 
obtained from Eq. (^) with 7 and Q.' replaced by 7k and Q' k (long dashed line) and with 7t and O' replaced by 7k and Q' k 
(short dash- long dashed line). 



Figure 12: Different expressions for the relative half-thickness h/r, monitored for a typical type II solution (a = 0, 
a — 0.5, M — 0.017, Lo = 0.9). The solid line describes the full expression (|42|). Dashed lines correspond to h/r as derived 
from Eq. (|43|), with 7 instead of 7k (short dashed) and 7t instead of 7k (long dashed). 



Figure 13: Parameter space showing contours of constant critical sound speed c s . The critical angular velocity Qf 1 
(bottom) and the corresponding critical specific angular momentum A* (top) are plotted for non-rotating (a = 0; left 
diagrams) and for a rotating (a — 0.95) black hole. Contours correspond to (from left to right) = 0.3, 0.15, 0.08, 0.04, 0.02, 
0.0. The c 2 s = 0.0 curve is identical with the Keplerian distributions Ak, Ok- In the S7-diagrams, short-dashed lines represent 
the limiting angular frequencies f2± for timelike trajectories, and long-dashed lines represent the angular velocity to of the 
LNRF. 



Figure 14: Variation of critical specific angular momenta A^ with the location of r s for a set of parameters E s . The 
upper diagram is for a = 0, the lower for a = 0.95. The outer dashed curves are the Keplerian distributions. Parameters E s 
are (from the top curves to the bottom) E s = 0.98, 1.00, 1.02, 1.03, 1.05 (dashed for a = 0), 1.10, 1.20, 1.50 respectively. 



Figure 15: Variation of the critical energy E s with the location of r s for a set of parameters A s . The lower dashed curves 
are the Keplerian distributions. For a — (upper diagram) parameters As are (from the top curve to the bottom) A = 0, 2, 3, 
3.13 (= A*; dashed), 3.4, 3.6, 3.8, 4 and 4.5 respectively. For a = 0.95 (lower diagram) parameters A s are (from the top curve 
to the bottom) A s = 0, 1.9, 2, 2.03 (= A*; dashed), 2.1, 2.2, 2.3, 2.4 and 2.5 respectively. Further plotted (dashed lines) are 
the curves of extrema -E max and -B m i n . 
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